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Abstract. Because of GAIA's estimated harvest of ~ 10 5 eclipsing bi- 
,-0 ■ naries (Munari et al. 2001) automative procedures for extracting physical 

P^P | parameters from observations must be introduced. We present prelimi- 

nary results of automative scripting applied to 5 eclipsing binaries, for 
which photometric and radial velocity observations were taken from liter- 
ature. Although the results are encouraging, extensive testing on a wider 
', sample has to be performed. 
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1. Introduction 



Q ' Eclipsing binaries have solid grounds in today's astrophysics; due to relatively 

simple geometry modeling allows astronomers to extract physical parameters 
such as masses and radii from observed light curves and radial velocity curves 
in absolute units, which is very difficult (if at all possible) for other kinds of 
Q-f celestial objects. Modeling approaches are developing rapidly with increasing 

q ■ computing power of PCs. Most widely used modeling algorithms are WD (Wil- 

son & Devinney 1971), WINK (Wood 1971), FOTEL (Hadrava 1986) and others. 
However, what most of these models lack is the level of automation required for 
scanning missions. GAIA will retrieve data for ~ 10 5 eclipsing binaries to 15 m 
(ESA-SCI(2000)4) and fully automative approaches are absolutely necessary. 

■ 2. Automative approach 

The fact that most of the existing modeling software lacks automation isn't a 
simple oversight, but is a consequence of few inevitable problems one faces while 
seeking a solution (Wilson 1998): 1. Extraction of physical parameters of eclips- 
ing binaries is a highly non- linear inverse problem and is as such subjected to 
modeling degeneracies. 2. Solution convergence may be slow or even fail because 
of shallow, wide minima in parameter space, numerical algorithms' inefficiency 
and non-linearity along with strong parameters' correlations. 3. The choice of 
parameters to be adjusted changes from iteration to iteration based on common 
sense and it is difficult to assess in advance how a certain change in parameter 
values changes the credibility of the entire solution. 

Bearing this in mind we use Wilson-Devinney model for basis and propose 
an automative approach, which is subject to the following main principles: (a) 
The underlaying formalism of solution-seeking shouldn't be changed in any way. 
(b) There is more to astrophysics than parameter estimation, so any automative 
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processing should be limited only to obtain coarse solutions, (c) The successful- 
ness should be tested on wide and diverse types of eclipsing binaries, (d) The 
time cost of modeling analysis should be considerably reduced. 

We refer to this process as automative scripting. Based on photometric 
and spectroscopic observations, with the assumption that the orbital period P, 
the epoch HJDo and the type of an eclipsing binary (detached, semi-detached, 
contact or over-contact) are known, we follow the predefined pattern of what 
modeling actions to perform in a particular order. Such a pattern is a script. 
Our goal is to create such a script (or a set of scripts) , that would lead us to a 
coarse modeling solution without any interactive approach. 

3. Preliminary test results 

We have performed some preliminary testing of automative scripting on dozen 
different eclipsing binaries taken from literature. Fig. 1 shows the results of 
automative scripting performed on 5 stars given in Table 1. The measurements 
consist of 100-1600 data points per filter in photometry and 10-100 data points 
in radial velocities. A typical time scale for the complete process of automative 
scripting on a 1.4GHz processor is ~35 seconds. Table 2 contains basic extracted 
parameters, which may be compared against adopted values given in Table 3. 



Star Name 


Type 


Available data 


References 


BH Vir 


Detached 


B, V, u, b, v, y, 


Clement et al. (1997) 






both RVs 




GK Dra 


Detached, 


B, V, both RVs 


Zwitter (2002) 




<5-Sct 2 ndar y 




TY Boo 


Over-contact 


B,V, both RVs 


Milone et al. (1991) 


UV Leo 


Detached 


B, V, both RVs 


Frederic & Etzel (1996) 


V505 Per 


Detached 


b T ,H p ,v T , both RVs 


Munari et al. (2001) 



Table 1. A list of chosen test stars that were used to assess the 
automative scripting successfulness. Notation remarks: U, B and V 
stand for Johnson, u, b, v, y for Stromgren, br, vt and H p for Ty- 
cho/Hipparcos and Ic for Cousins' filters. RVs are radial velocities. 

Few remarks should be made: 1. The solutions were reached without any 
human intervention. 2. The solution degeneracy is reduced by using the method 
of multiple subsets (Wilson 1998) instead of exclusive use of differential correc- 
tions. 3. Further degeneracy elimination is possible because of the presence 
of radial velocities, which along with photometric data uniquely determine the 
semi-major axis and both components' masses. Because of GAIA's RVS capa- 
bilities we expect tremendous improvements to the field of eclipsing binaries. 

4. Discussion 

Note the discrepancy between calculated and adopted potentials fii and 0,2 
in Tables 2 and 3: this is a consequence of the degeneracy between the in- 
clination and stellar radii (oc $7 _1 ), the direct example of how blindly trusting 




Figure 1. Light curves (left) and radial velocity curves (right) of the 
5 stars given in Table 1: BH Vir, GK Dra, TY Boo, UV Leo and 
V505 Per, respectively from top. The independent variable is in all 
cases phase and the dependent variable is relative flux for light curves 
and normalized radial velocities to 2na/P for radial velocity curves. 
Measurements are marked with filled circles and the synthetic solution 
from Table 2 with a solid line. 
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automatically retrieved solutions may be dangerous. Since we can recognize 
this problem, we may avoid it by using more adequate scripts, but there may be 
many caveats we must still anticipate. 



Star Name 


a[R Q ] 


q 


T 2 /Ti 




Oi 


n 2 


w 7 [km/s] 


BH Vir 


4.7(3) 


0.90(0) 


0.92(0) 


85.9(1) 


5.2(0) 


4.7(0) 


-0.2(8) 


GK Dra 


28.9(4) 


1.26(0) 


1.01(0) 


83.9(2) 


10.5(1) 


10.5(1) 


1.7(12) 


TY Boo 


2.0(2) 


2.15(1) 


1.00(0) 


76.4(1) 


4.0(4) 


4.0(4) 


-43.8(15) 


UV Leo 


4.0(3) 


1.00(1) 


0.97(1) 


84.1(1) 


5.1(0) 


4.4(0) 


12.8(6) 


V505 Per 


15.1(1) 


0.98(1) 


0.99(0) 


87.0(2) 


11.5(2) 


11.2(2) 


-0.7(6) 



Table 2. Computed values of basic stellar parameters using automa- 
tive scripting. Standard deviation is given in parentheses as a measure 
of accuracy of the last decimal place. These are modeling standard 
deviations and may be larger in physical sense (Wilson 1998). 



Star Name 


a[R Q ] 


q 


T 2 /T 1 




Oi 


^2 


v-y[km/s] 


BH Vir 


4.7 


0.89 


0.91 


87.5 


4.7 


4.8 


0.0 


GK Dra 


29.0 


1.26 


0.98 


85.8 


14.1 


12.2 


3.8 


TY Boo 


2.32 


2.13 


1.00 


77.5 


5.4 


5.4 


-38.7 


UV Leo 


3.91 


0.96 


0.98 


82.6 


4.6 


4.3 


13.0 


V505 Per 


15.0 


0.98 


0.98 


87.8 


11.5 


11.2 


0.0 



Table 3. Adopted values of basic stellar parameters from references 
given in Table 1. 

The results demonstrated in the previous section are encouraging; how- 
ever, one has to keep a fair amount of scientific scepticism about the adopted 
automation philosophy. This study was performed with an extremely undersam- 
pled statistics. A thorough study over large numbers of diverse eclipsing binaries 
is compulsory to gain confidence in the automation procedure. At this time we 
are still uncertain about the universality of the developed algorithm. We may 
be facing systematic errors due to conformity of our script: it may invoke the 
existence of a preferred subspace within the parameter space, a dire consequence 
that would compromise all retrieved solutions. 
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